Search results

The use of the prominent FDTD method for the time domain solution of electromagnetic wave propagation past devices with small geometrical details can require very fine grids and can lead to very important computational time and storage. The purpose is to develop a numerical method able to handle possibly non‐conforming locally refined grids, based on portions of Cartesian grids in order to use existing pre‐ and post‐processing tools.

A Discontinuous Galerkin method is built based on bricks and its stability, accuracy and efficiency are proved.

It is found to be possible to conserve exactly the electromagnetic energy and weakly preserves the divergence of the fields (on conforming grids). For non‐conforming grids, the local sets of basis functions are enriched at subgrid interfaces in order to get rid of possible spurious wave reflections.

Although the dispersion analysis is incomplete, the numerical results are really encouraging it is shown the proposed numerical method makes it possible to handle devices with extremely small details. Further investigations are possible with different, higher‐order discontinuous finite elements.

This paper can be of great value for people wanting to migrate from FDTD methods to more up to date time‐domain methods, while conserving existing pre‐ and post‐processing tools.

The use of the prominent finite difference time‐domain (FDTD) method for the time‐domain solution of electromagnetic wave propagation past devices with small geometrical details can require very fine grids and can lead to unmanageable computational time and storage. The purpose of this paper is to extend the analysis of a discontinuous Galerkin time‐domain (DGTD) method (able to handle possibly non‐conforming locally refined grids, based on portions of Cartesian grids) and investigate the use of perfectly matched layer regions and the coupling with a fictitious domain approach. The use of a DGTD method with a locally refined, non‐conforming mesh can help focusing on these small details. In this paper, the adaptation to the DGTD method of the fictitious domain approach initially developed for the FDTD is considered, in order to avoid the use of a volume mesh fitting the geometry near the details.

Based on a DGTD method, a fictitious domain approach is developed to deal with complex and small geometrical details.

The fictitious domain approach is a very interesting complement to the FDTD method, since it makes it possible to handle complex geometries. However, the fictitious domain approach requires small volume elements, thus making the use of the FDTD on wide, regular, fine grids often unmanageable. The DGTD method has the ability to handle easily locally refined grids and the paper shows it can be coupled to a fictitious domain approach.

Although the stability and dispersion analysis of the DGTD method is complete, the theoretical analysis of the fictitious domain approach in the DGTD context is not. It is a subject of further investigation (which could provide important insights for potential improvements).

This is believed to be the first time a DGTD method is coupled with a fictitious domain approach.

The purpose of this paper is to develop a time implicit discontinuous Galerkin method for the simulation of two‐dimensional time‐domain electromagnetic wave propagation on non‐uniform triangular meshes.

The proposed method combines an arbitrary high‐order discontinuous Galerkin method for the discretization in space designed on triangular meshes, with a second‐order Cranck‐Nicolson scheme for time integration. At each time step, a multifrontal sparse LU method is used for solving the linear system resulting from the discretization of the TE Maxwell equations.

Despite the computational overhead of the solution of a linear system at each time step, the resulting implicit discontinuous Galerkin time‐domain method allows for a noticeable reduction of the computing time as compared to its explicit counterpart based on a leap‐frog time integration scheme.

The proposed method is useful if the underlying mesh is non‐uniform or locally refined such as when dealing with complex geometric features or with heterogeneous propagation media.

The paper is a first step towards the development of an efficient discontinuous Galerkin method for the simulation of three‐dimensional time‐domain electromagnetic wave propagation on non‐uniform tetrahedral meshes. It yields first insights of the capabilities of implicit time stepping through a detailed numerical assessment of accuracy properties and computational performances.

In the field of high‐frequency computational electromagnetism, the use of implicit time stepping has so far been limited to Cartesian meshes in conjunction with the finite difference time‐domain (FDTD) method (e.g. the alternating direction implicit FDTD method). The paper is the first attempt to combine implicit time stepping with a discontinuous Galerkin discretization method designed on simplex meshes.

The purpose of this paper is to analyse algorithms for fluid‐structure interaction (FSI) from a purely algorithmic point of view.

First of all a 1D model problem is selected, for which both the fluid and structural behavior are represented through a minimum number of parameters. Different coupling algorithm and time integration schemes are then applied to the simplified model problem and their properties are discussed depending on the values assumed by the parameters. Both exact and approximate time integration schemes are considered in the same framework so to allow an assessment of the different sources of error.

The properties of staggered coupling schemes are confirmed. An insight on the convergence behavior of iterative coupling schemes is provided. A technique to improve such convergence is then discussed.

All the results are proved for a given family of time integration schemes. The technique proposed can be applied to other families of time integration techniques, but some of the analytical results need to be reworked under this assumption.

The problems that are commonly encountered in FSI can be justified by simple arguments. It can also be shown that the limit at which trivial iterative schemes experience convergence difficulties is very close to that at which staggered schemes become unstable.

All the results shown are based on simple mathematics. The problems are presented so to be independent of the particular choice for the solution of the fluid flow.

This work is concerned with the development and the numerical investigation of a hybridizable discontinuous Galerkin (HDG) method for the simulation of two‐dimensional time‐harmonic electromagnetic wave propagation problems.

The proposed HDG method for the discretization of the two‐dimensional transverse magnetic Maxwell equations relies on an arbitrary high order nodal interpolation of the electromagnetic field components and is formulated on triangular meshes. In the HDG method, an additional hybrid variable is introduced on the faces of the elements, with which the element‐wise (local) solutions can be defined. A so‐called conservativity condition is imposed on the numerical flux, which can be defined in terms of the hybrid variable, at the interface between neighbouring elements. The linear system of equations for the unknowns associated with the hybrid variable is solved here using a multifrontal sparse LU method. The formulation is given, and the relationship between the considered HDG method and a standard upwind flux‐based DG method is also examined.

The approximate solutions for both electric and magnetic fields converge with the optimal order of p+1 in L₂ norm, when the interpolation order on every element and every interface is p and the sought solution is sufficiently regular. The presented numerical results show the effectiveness of the proposed HDG method, especially when compared with a classical upwind flux‐based DG method.

The work described here is a demonstration of the viability of a HDG formulation for solving the time‐harmonic Maxwell equations through a detailed numerical assessment of accuracy properties and computational performances.

The simulation of the fluid–solid interaction (FSI) problem is important for both academic studies and engineering applications. However, the numerical approach for simulating the FSI problems is a great challenge owing to the large discrepancy of material properties and inconsistent description of grid motion between the fluid and solid domains. The difficulties will be further increased if there are multiple materials in the fluid region. In these complicated applications, interface reconstruction, multi-material advection and FSI must be all taken into account. This paper aims to present an effective integrated work of multi-material arbitrary Lagrangian Eulerian (MMALE) method, finite element (FE) method and the continuum analogy method to simulate the complex FSI problems involving multi-material flow. The coupled method is used to simulate the three-dimensional CONT test and the blast-plate interaction. The numerical results show good agreement with the benchmark and the experiment data, which indicates that the presented method is effective for solving the complicated FSI problems.

MMALE and FE methods are used to simulate fluid and solid regions, respectively. The interfacial nodes of fluid and solid are required to be coincident in the whole simulation so the interacted force can be easily and accurately calculated. To this end, the continuum analogy method is used in the rezoning phase.

The coupled method is used to simulate the three-dimensional CONT test and the blast-plate interaction. The numerical results show good agreement with the benchmark and the experiment data, which indicates that the presented method is effective for solving the complicated FSI problems.

To the best of the authors’ knowledge, this is the first time that the ALE method, moment of fluid interface reconstruction method, continuum analogy method and the FE method are combined to solve complicated practical problems.

This study contributes to the literature on local administrations' debt and attempts to answer the following research questions: (1) What effects do swaps produce on regions' debt? (2) Have swaps been used to finance discretionary debt?

The paper investigates the debt burden as influenced by economic, financial and political variables and forces with panel data techniques, and tests whether swaps have been used to financing debt due to unfunded expenditures.

Panel data results of 15 Italian regions over the 2007–2014 period shows that regions with higher debt exhibited a higher interest rate exposure and have employed derivatives hoping to counterbalance the reduced resources received from the central state, in line with other European countries' experience (i.e. France and Greece).

The scarcity of official data and information on swaps has limited the empirical investigations in the literature but did not reduce the losses of local administrations.

This study creates the first database on swaps purchased by Italian regions to investigate their impact on their debt. Results show that highly indebted regions with reduced funds from the central state and diminished local resources are more likely to use swaps to fund their debt. Italian regions heavily depended on long-term debt to finance their non-healthcare services, rather than current revenues; swaps have been used to finance discretionary (non-healthcare) debt.

This paper aims to present a fluid-structure coupling partitioned scheme involving rigid bodies supported by spring-damper systems. This scheme can be used with already existing fluid flow solvers without the need to modify them.

The scheme is based on a modified Broyden method. It solves the equations of solid body motion in which the external forces coming from the flow are provided by a segregated flow solver used as a black box. The whole scheme is implicit.

The proposed partitioned method is stable even in the ultimate case of very strong fluid–solid interactions involving a massless cylinder oscillating with no structural damping. The overhead associated with the coupling scheme represents an execution time increase by a factor of about 2 to 5, depending on the context. The scheme also has the advantage of being able to incorporate turbulence modeling directly through the flow solver. It has been tested successfully with URANS simulations without wall law, thus involving thin high aspect-ratio cells near the wall.

Such problems are known to be very difficult to solve and previous studies usually rely on monolithic approaches. To the authors' knowledge, this is the first time a partitioned scheme is used to solve fluid–solid interactions involving massless components.

The purpose of this study is twofold: first, to derive a consistent model free-surface runup flow problems over deformable beds. The authors couple the nonlinear one-dimensional shallow water equations, including friction terms for the water free-surface and the two-dimensional second-order solid elastostatic equations for the bed deformation. Second, to develop a robust hybrid finite element/finite volume method for solving free-surface runup flow problems over deformable beds. The authors combine the finite volume for free-surface flows and the finite element method for bed elasticity.

The authors propose a new model for wave runup by static deformation on seabeds. The model consists of the depth-averaged shallow water system for the water free-surface coupled to the second-order elastostatic formulation for the bed deformation. At the interface between the water flow and the seabed, transfer conditions are implemented. Here, hydrostatic pressure and friction forces are considered for the elastostatic equations, whereas bathymetric forces are accounted for in the shallow water equations. As numerical solvers, the authors propose a well-balanced finite volume method for the flow system and a stabilized finite element method for elastostatics.

The developed coupled depth-averaged shallow water system and second-order solid elastostatic system is well suited for modeling wave runup by deformation on seabeds. The derived coupling conditions at the interface between the water flow and the bed topography resolve well the condition transfer between the two systems. The proposed hybrid finite volume element method is accurate and efficient for this class of models. The novel technique used for wet/dry treatment accurately captures the moving fronts in the computational domain without generating nonphysical oscillations. The presented numerical results demonstrate the high performance of the proposed methods.

Enhancing modeling and computations for wave runup problems is at an early stage in the literature, and it is a new and exciting area of research. To the best of our knowledge, solving wave runup problems by static deformation on seabeds using a hybrid finite volume element method is presented for the first time. The results of this research study, and the research methodologies, will have an important influence on a range of other scientists carrying out research in related fields.

The purpose of this paper is to present a time domain discontinuous Galerkin (DG) approach for modeling wideband frequency dependent surface impedance boundary conditions.

The paper solves the Maxwellian initial value problem in a computational domain, which is spatially discretized by the higher order DG method. On the boundary of the computational domain the paper applies a suitable impedance boundary condition (IBC). The frequency dependency of the impedance function is modeled by auxiliary differential equations (ADE).

The authors will study the resonance frequency and the Q factor of different types of cavity resonators including lossy materials. The lossy materials are modeled by means of IBCs. The authors will compare the results with analytical results, as well as numerical results obtained by direct calculations where lossy materials are included explicitly into the numerical model. Several convergence studies are performed which demonstrate the accuracy of the approach.

Modeling of frequency dependent boundary conditions in time domain with finite difference time domain method (FDTD) method is considered in numerous papers, as well as in frequency domain finite element method (FEM), and in a few papers also time domain FEM. However, FDTD method is only first order accurate and fails in modeling of complicated surfaces. FEM allows for high order accuracy, but time domain modeling is numerically extremely expensive. In frequency domain, broadband modeling of frequency dependent boundary conditions requires several simulations as opposed to the time domain, where a single simulation is needed. The time domain DG method proposed in this paper allows to overcome the difficulties. The authors introduce a broadband surface impedance formulation based on the ADE approach for the higher order DG method.